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ANALYSIS OF A BORON-CARBIDE-DRUM-CONTROLLED 
CRITICAL REACTOR EXPERIMENT 
by Wendell Mayo 
Lewis Research Center 

SUMMARY 

Fast-spectrum nuclear reactors as heat sources for the generation of electric power 
in space are part of the technology program at the Lewis Research Center. A neutron 
absorber, such as boron carbide enriched in boron- 10, incorporated in control drums is 
one type of reactivity control method that is being considered. In order to validate meth- 
ods and cross sections used in the neutronic design of these reactor concepts, critical 
experiments were performed at Atomics International to check design calculations. The 
experiments incorporated materials for high -temperature, long-life reactors such as 
uranium -235 nitride fuel, lithium-7 coolant, and refractory -metal structure and reflec- 
tors. The critical assembly diameter is close to the diameter considered for the power 
reactor concepts. The fuel -pin-in-honeycomb-lattice core design and molybdenum radial 
and axial reflectors were also simulated in the critical experiments. The reactor was 
controlled by rotating control drums containing boron carbide. 

The multigroup two-dimensional transport analysis gave generally satisfactory re- 
sults. The calculated multiplication factor for the most detailed spatial calculation was 
only 0. 7 -percent Ak too high. The calculated reactivity worth of the control drums was 
$11. 61 compared to $11. 58 measured by the inverse kinetics method and $11. 98 meas- 
ured by the inverse counting method. Calculated radial and axial power distributions 
were in good agreement with experiment. 


INTRODUCTION 

Fast-spectrum reactor concepts that use various reactivity control methods have 
been considered in the technology program at the Lewis Research Center. The nuclear 
reactors would serve as heat sources for systems generating electric power in space. 
Reversible hydriding of yttrium zones (ref. 1), axially translating radial reflectors 



(ref. 2), rotating poisoned control drums (ref. 3), and rotating fueled control drums 
(refs. 4 and 5) are examples of control methods previously considered. Performance 
goals established were generally for about 2. 2 thermal megawatts for 50 000 hours with 
a lithium -7 coolant outlet temperature of 1222 K. 

The conclusions in reference 3 indicate that reactivity control using boron- 10 carbide 
control drums located outside the pressure vessel might be marginal. Since that report 
was written (1969), several other design concepts using boron-10 carbide for control have 
been examined. These include the use of in-core control rods, the use of control rods in 
conjunction with control drums, and the use of control drums located in dry wells inside 
the pressure vessel. In each case, the major consideration was to provide adequate re- 
activity control for the long-life (50 000 hr), 2. 2-thermal megawatt reactor requirements. 

More recently, as part of an experimental program (ref. 6) under contract NAS 
3-12982 and an extension of that program (ref. 7) under contract NAS 3-14421, a critical 
assembly using boron carbide control drums was built and operated by Atomics Inter- 
national. Preliminary calculations of these experiments indicated an underestimation of 
the measured control -drum worth. With this discrepancy in mind, a reexamination of the 
methods and cross sections used in the calculation of the critical assembly and in the pre- 
vious design concept calculations of boron-carbide -controlled reactors (such as in ref. 3) 
is in order. 

Although a comprehensive series of experiments were performed at Atomics Inter- 
national, the present work is limited to those associated with the boron -carbide -drum - 
controlled reactor. This report presents an analysis of the reactor critical experiments 
to determine possible sources of the discrepancy in control -drum reactivity calculations 
mentioned previously. 


DESCRIPTION OF EXPERIMENTS 

Complete details of the critical assembly design and experiments are given in ref- 
erences 6 and 7 so that only a brief description pertinent to this report is given here. 

Figure 1 shows the fuel element used in the critical assembly. The fuel element has 
molybdenum end reflectors. The reflector piece and the ceramic-like lithium-7 nitride 
are eccentrically sealed between the honeycomb tube and the fuel tube. The cylindrical 
molybdenum reflector piece can be taken out of the fuel tube after removal of the 
jackscrew - neoprene -seal assembly, thus providing access to the fuel region for fuel 
loading adjustments. Uranium metal rods (93. 145-percent U) are used for the fuel. 
Three sizes are used so that relatively fine adjustments can be made to assure fuel load- 
ing uniformity between fuel elements. The fuel rod bundle is centered inside the fuel tube 
by tantalum centering rings at each end of the fuel bundle. The 0. 051 -centimeter offset 
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between the fuel tube and the honeycomb tube is for the purpose of measuring radial 
compression-expansion effects resulting from fuel displacement, as would occur with 
temperature changes in the design core. The normal (null) position is such that the mark 
on the end cap (see fig. 1) is perpendicular to the radius vector to the centerline of the 
core. 

The fuel elements are positioned in a lattice (with a pitch of 2. 215 cm), as shown in 
figure 2, by grid plates with holes into which the fuel -element end fittings are inserted. 
The lattice forms a six -pointed star shape composed of 181 fuel elements. Between the 
star -points are located the control drums. As shown in figure 2 the boron carbide is in- 
side aluminum cannisters bolted to the massive molybdenum drum center piece. Oppo- 
site the boron carbide is another aluminum cannister filled with molybdenum plates. The 
stationary and movable radial reflector pieces are also of molybdenum. The design is 
such that four of the six radial reflector pieces can fall away from the core for scram- 
ming the reactor. The tantalum (simulating a pressure vessel) attached to the radial re- 
flector pieces is therefore also segmented. 

Measurements included the critical loading (or excess reactivity near critical), the 
radial and axial power distributions, and the worth of single control drums and of all six 
control drums. Details of the measurements and procedures used are fully described in 
references 6 and 7. 


ANALYSIS PROCEDURE 

In this section the geometric models and procedures used in the analysis of the ex- 
periments are described. Sufficient detail is provided so that the same calculations can 
be repeated, using updated cross sections as they become available, if desired. 

Three separate energy-group structures are used, as shown in table I. Previous 
calculations have been done primarily with the four -energy -group cross sections, with 
occasional checks made by using the 13 -group set. Generally, the four -group set has 
been found adequate for preliminary calculations and for surveys. A 26 -group set was 
employed for the critical experiment analysis to check the adequacy of the fewer -group 
treatment. 

Two sets of cross sections were used. One set, designated GAM, is fairly old and 
was used in nearly all previous calculations during the past few years. The second set, 
labeled GAM(ENDF), incorporates a number of isotopes from the ENDF/B cross-section 
library as processed by the SUPERTOG program (ref. 8) for use with the GAM -II pro- 
gram (ref. 9). This program was used for spectrum averaging all cross sections used in 
this study. Table II lists the ENDF/B material numbers used in GAM -II. With the 
GAM(ENDF) set, the older GAM isotopes were used to supplement the list in table II as 
required. 
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Figure 3(a) shows one of the 1/4 -core calculation models used in the XY calcula- 
tions. The XY grid was superimposed over the experimental geometry shown in fig- 
ure 2 in such a manner as to represent the important regions. A certain amount of ho- 
mogenizing is done to reduce mesh interval requirements. Briefly, the calculation model 
geometry is as follows: 

(1) The fuel region is based on cells. The materials are smeared within the cell. 

(2) The aluminum cans around the boron carbide powder and the molybdenum plates 
are smeared into the reflector region, which includes all materials inside the pressure 
vessel except the boron carbide regions and the fuel regions. 

(3) In order to conserve mesh intervals, the exact area of the pressure vessel is not 
preserved; however, the atom density of the tantalum is adjusted accordingly to preserve 
the correct number of atoms. 

(4) For one set of calculations, which incorporate a void explicitly between the con- 
trol drums and the core fuel cells, the reflector materials and area are adjusted some- 
what, though the total number of atoms of the reflector materials is preserved (fig. 3(b)). 

The input axial and radial mesh specifications and the material numbers by region as 
used in the XY spatial calculations are listed in the appendix just as used in the DOT- 
IIW programs (ref. 10). Also listed are the atom densities for each region (table VII). 
The weights and dimensions of each piece of the critical assembly as described in refer- 
ences 6 and 7 were used in determining the atom densities and regions. 

In order to account for the end reflector and axial effects in the calculations, an RZ 
model was also established. Figure 4 shows the RZ model used; and table m lists the 
volume percents of constituents in each region, along with consistent densities of the ma- 
terials, so that the atom densities used in the calculations can be reproduced. In the ra- 
dial direction, at Z = 0, the RZ model is consistent with the XY model from which it 
was derived. 

Information obtained from the RZ calculations consists of axial power distributions 
and the effective height for use in determining transverse leakage ’’cross sections” for 
use in the XY calculations. The transverse leakage cross section, is defined as 



where Z^ r is the macroscopic transport cross section for group g and region i, and 
g 

H eff is an ef f ective height (core length + reflector savings) determined by running one- 
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dimensional radial calculations (Z = 0) and varying the height until the same multiplica- 
tion factor as given by the RZ calculation is obtained. For a given cross-section set 
and quadrature order, is assumed to be constant whether the poisoned portion of 
the control drums are rotated fully in or fully out. Previous calculations have indicated 
this is a good assumption. The XY calculations are used to obtain the radial power dis- 
tributions, the effective multiplication factor, and the reactivity control worth of the 
drums. 


COMPARISON OF EXPERIMENTS AND CALCULATIONS 

The fuel loading in the critical assembly was 143. 797 kilograms (corrected for the 
protective coating) of 93. 145 -percent enriched uranium. This fuel loading resulted in an 
excess reactivity of $1. 606 with the poisoned portion of the control drums rotated fully 
out. With an effective delayed neutron fraction of 0. 0067, this excess reactivity corre- 
sponds to a multiplication factor k of 1. 011; the calculations in this section are to be 
compared with this value. 


Effect of Low-Order to High-Order Calculations 

Table IV lists calculations made to determine the effect of varying the number of en- 
ergy groups, the quadrature order, and the order of the Legendre expansion of the scat- 
tering cross sections. GAM cross sections were used, and both the RZ geometry and 
XY geometry results are tabulated for the "poison-out" case. The XY geometry cor- 
responds to the second set in the appendix with 36X by 34Y intervals. The trend in mul- 
tiplication factor generally follows that observed previously (refs. 5, 11, and 12) in that 
the higher order calculations give lower multiplication factors. Little benefit is ob- 
served in going beyond the 13-group S^P^ treatment. The results have been rounded to 
the third decimal, corresponding to about ±0. 05-percent-Ak/k (±7. 5 cents) accuracy. 
Because of slight variations in the calculated effective height, the XY results show less 
change between low -order and high -order calculations. The XY results are the ones to 
compare with the experimental multiplication factor of 1.011 since the XY calculations 
incorporate the axial effects (through the effective height) and account for the azimuthal 
asymmetry more accurately than the RZ model can. Thus the higher order calculations 
give a multiplication factor 1. 2 -percent Ak too large when this geometric model is used. 
Another XY model, which more closely matches the experiment in that an explicit void 
region appears between the control drums and the core, was also used (fig. 3(b), geom- 
etry 4 in the appendix). This more-detailed model gave better agreement; k was only 
0.7 -percent Ak too large. 
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Effect of Using GAM and GAM(ENDF) Cross Sections 

The calculations were made in the manner described in the preceding section, except 
that GAM(ENDF) cross sections (table II) were substituted for the GAM set. The GAM- 
(ENDF) set gave a multiplication factor 0.9 -percent Ak larger than the GAM set for the 
13 -group calculation procedure. The 26 -group set showed a similar tendency to 

calculate larger multiplication factors. Since uranium -235 is an important and major 
constituent in the critical experiments, an independent comparison was made between the 
GAM and GAM(ENDF) cross sections for the fuel. 

Reference 13 describes a GODIVA-like spherical critical experimental performed at 
the Oak Ridge National Laboratory. Calculation of the experiment using GAM cross sec- 
tions (SgPp 26 group) gave a multiplication factor 0. 6-percent Ak larger than the ex- 
perimental value. (This is fairly consistent with the 0. 7-percent Ak overestimate dis- 
cussed in the preceding section. ) Use of the GAM(ENDF) cross sections (material 1044 
235 

for U) indicated a multiplication factor 0. 8 -percent Ak larger than when the GAM set 
was used, which is again consistent with the results for the boron-carbide -drum - 
controlled experiment. Thus, we are led to suspect the GAM(ENDF) cross sections for 
uranium -235 (material 1044). The ENDF/B uranium cross sections have been updated 
but are not available as yet for use with the GAM program. It should be pointed out that, 
as in most calculations of the complexity of those in this study, off -setting discrepancies 
in cross sections can lead to good agreement with experimental results. In this study, 
the boron- 10 cross sections in GAM are clearly in error (discussed in the next section); 
the GAM(ENDF) set may also be wrong but to a lesser extent. Similarly, we can infer 
probable errors in the uranium -235 cross sections based on the calculation of these two 
critical experiments. However, an element -by -element substitution for the other mate- 
rials was not done, which precludes any judgement as to their accuracy. 


Drum Control Worth 

The reactivity worth of the control drum was measured in two different ways. One 
measurement used was the inverse kinetics method to measure the worth of each control 
drum. The average worth per drum was $1. 93. Assuming no adjacent drum interac- 
tions, the worth of six drums would be $11. 58. The inverse counting method was also 
used to measure the worth of all six drums ganged. This method is not considered highly 
accurate but gave a reactivity worth of $11. 98, which is in fair agreement with the in- 
verse kinetics measurement. Both methods are explained in detail in references 6 and 7. 

Initial calculations for the reactivity worth of the control drums were too low for two 
reasons: 
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(1) The old GAM cross sections for boron 10 had too little absorption in the energy 
range of importance - between about 100 keV and 2 MeV. This was suspected (ref. 14) 
and the ENDF/B cross-section library (material 1009) reflects significantly larger ab- 
sorption cross sections. Table V shows the GAM and GAM(ENDF) cross sections for the 
26 -energy -group structure. 

(2) Initial calculations were done with reflector material smeared in between the con- 
trol drums and the active core (XY geometries 1 and 2 in the appendix). This reflector 
material effectively shields the control drum and reduces control swing. The XY geom- 
etries 3 and 4 in the appendix incorporate an explicit void region between the drums and 
the core which more closely matches the experimental geometry. Table VI gives results 
of the various calculations. The control swing is not very sensitive to the order of the 
calculations; compare cases 1, 2, and 3 and compare cases 5 and 7. Also note that the 
use of the boron-10 cross sections from ENDF/B (cases 4, 5, and 7) have control swings 
about $1 larger than those calculated by using the older GAM cross sections for boron- 10 
(cases 1 to 3). Case 6 is special in that it was performed by using the XY geometry 
setup to show the explicit void region between control drums and core. The reactivity 
worth of $11.61 is in good agreement with the experimental values. 


Axial Power Distributions 

The relative axial power distribution calculated with the RZ model is shown in fig- 
ure 5, along with the experimental points (ref. 7). The measurements were made by 
gamma-scanning small -diameter fuel wires. The experimental points at three radial 
positions were averaged to obtain a single value at each axial location. Several experi- 
mental values (about 4 percent of the total measured points) that showed unexplained ir- 
regularities were discarded in the averaging process. Also, since the RZ calculation 
assumes symmetry at the axial midplane, the corresponding experimental points above 
and below the midplane were averaged to obtain the values plotted in figure 5. An inspec- 
tion of the data in reference 7 indicates that the separability of radial and axial flux solu- 
tions is a good assumption and that axial asymmetry was not significant. Thus, calcula- 
tion and experiment (fig. 5) are generally in good agreement. 


Radial Power Distributions 

Figure 6 shows a 30° sector of the critical assembly. The relative power density 
values for calculation and experiment are shown with the corresponding fuel elements and 
at two locations in the boron carbide region. Three axial points for each radial position 
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were measured and averaged to obtain the values shown in the figure. Again, a small 
number of unexplained irregular points were discarded in the averaging procedure. 
Agreement is fairly good for all points shown in figure 6. The 9 percent difference in the 
outer boron carbide region is the largest difference noted, but it is in a relatively unim- 
portant region; the multiplication factor would not be sensitive to this discrepancy. 


SUMMARY OF RESULTS 

A boron-carbide -drum -controlled fast-spectrum critical assembly was analyzed to 
check calculation methods and cross sections used in nuclear reactor design. The crit- 
ical assembly used enriched uranium fuel, tantalum structure, lithium and nitrogen in 
the form of lithium nitride, molybdenum reflectors, and boron carbide control drums; 
all these materials are being considered for possible use in space power reactor 
concepts. 

It was found that calculations using 13 energy groups were adequate for the 

analysis. Both GAM and GAM(ENDF) cross sections were used. In computing the mul- 
tiplication factor the GAM set was 0. 7-percent Ak high, while the GAM(ENDF) set was 
1. 6-percent Ak too large based on the most detailed two-dimensional geometric repre- 
sentation. An independent check of a small fast -spectrum critical reactor indicated that 
the larger error with the GAM(ENDF) cross sections was due to the GAM(ENDF) uranium - 
235 (material 1044). 

In calculating the reactivity worth of the control drums, it was found that the newer 
boron- 10 cross sections (ENDF material 1009) were preferable. When the GAM(ENDF) 
boron- 10 cross sections were used, the calculated reactivity worth of the control drums 
was $11. 61 compared to $11. 58 measured by the inverse kinetics method and $11.98 
measured by the inverse counting method. With the older GAM boron-10 cross sections, 
an underestimate of about $1 was observed. 

The importance of a detailed geometric model was confirmed in the calculation of the 
control worth of the drums. Explicitly representing the void region between drums and 
core increased the control swing by about $1 compared with a less -detailed model in 
which reflector material was smeared into that region. 

Radial and axial power distributions were in good agreement with the measurements. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 5, 1972, 

503-25. 
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APPENDIX - MATERIAL AND GEOMETRY SPECIFICATIONS FOR CALCULATIONS 


FOR B 4 C-DRUM-CONTROLLED REACTOR EXPERIMENT 

This appendix records the atom densities (table VII) for each material used in the 
XY calculations and provides card images of the geometric specifications as used by the 
DOT-IIW program for each of the four geometries. The geometries are referred to in 
the text by numerical order in which they appear in the appendix. Briefly, the designa- 
tions are 

(1) 37X by 35Y intervals, control -drum poison inserted, six materials (regions) 

(2) 34X by 36Y intervals, control -drum poison rotated full out, six materials 

(regions) 

(3) 37X by 35Y intervals, control -drum poison inserted, seven materials (regions), 

explicit void region between control drums and core 

(4) 44X by 45Y intervals, control -drum poison rotated full out, seven materials 

(regions), explicit void region between control drums and core 
The card images are as follows: 
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3 

4R 

2 


9R 

6 








22R 

2 

3R 

5 9R 

6 





22R 

2 


51 1R 

6 





20R 

2 

3R 

51 1R 

6 





20R 

2 


513R 

6 





17R 

2 

4R 

513R 

6 





17R 

2 


516R 

6 





6R 

4 

9R 

2 3R 

516R 

6 




6R 

4 

9R 

2 

518R 

6 




6R 

4 

7R 

2 3R 

518R 

6 




4R 

4 

7R 

2 3R 

520R 

6 




4R 

4 

7R 

2 

522R 

6 




3R 

4 

6R 

2 3R 

522R 

6 




3R 

4 

6R 

2 

524R 

6 




2R 

4 

5R 

2 3R 

524R 

6 




3R 

2 

5R 

526R 

6 





4R 

530R 

6 







37X BY 3 5 Y GEOMETRY - MESH AND MATERIAL SPECIFICATIONS - POISON IN 
WITH EXPLICIT VOIDS BETWEEN DRUMS AND CORE 

4* 0.0 


1.08 

2 . C 8 

3.3 


4.20 

4.85 

5.40 

6 » 20 

6.7 

7.7 


8.7 

9.7 

10.6894 

li.O 

1 1 . 3 

11.6712 


12.5 

13.41 

14.20 

14.4 

15.60 

16.3 


17.6 

18.25 

19.20 

19.80 

20.1 

21.62 


23.00 

23.7 

24.80 

<:5.S0 

26.50 

27.20 


27.80 

28.20 

28.6 

29.291 







2* 

o 

• 

o 






1 .08 

2.215 

3.3 


3.45 

4.4 

4.80 

5.54 

5.89 

6.4 


7.20 

7.5 

8. 80 

9.50 

10.51 

12.2 


13.25 

13.8 

14. 3 

14.6 

15.4 

15.7 


16.5 

17.6 

18.812 

19.1844 

21.70 

22.40 


23.60 

24.758 

25.90 

^'3.90 

2Y.c>0 

28.20 


28.6 

29.291 


8$ 







26R 

1 10R 2 


5 




23R 

113R 2 


5 . 




2 1 R 

l i 5R 2 


5 





Ui O' O' 



19R 

.1 3R 

714R 

2 

5 


19R 

1 

7 5R 

31 1R 

2 

5 

1 6R 

1 3R 

7 5R 

312R 

2 

t> 

1 6R 

1 

7 7R 

312R 

2 

5 

12R 

1 3R 

7 

2 

7 7R 

312R 


5 





12R 

1 3R 

2 

7 7R 

311R 

2 3R 

12R 

1 3R 

2 8R 

31 1R 

2 

3 2R 

12R 

1 2R 

2 

7 7R 

312R 

2 

2R 

6 





12R 

I 

2 

7 8R 

312R 

2 

2R 

6 





12R 

1 

7 8R 

31 1R 

2 3R 

5 2R 

I2R 

1 

7 8R 

3 1 1R 

2 

5 4R 

12R 

1 

7 7R 

312R 

2 

5 4R 


ZilR 

1 2R 

7 5R 

313R 

2 

4R 

6 





2R 

2 

7 9R 

1 3R 

2 3R 

312R 

3R 

5 4R 

6 




2R 

4 

7 9R 

1 3R 

2 3R 

3 12R 


5 6R 

6 




3R 

4 2R 

7 7R 

1 4R 

2 

313R 


5 6R 

6 




4R 

4 

7 7R 

1 4R 

2 

313R 


5 6R 

6 




5R 

4 3R 

7 4R 

118R 

2 

5 6R 

6R 

4 

7 

2 4R 

1 16R 

2 3R 

6R 

6 





7R 

4 3R 

2 2R 

1I6R 

2 

5 8R 

8R 

420R 

2 

5 8R 

6 


8R 

4i7R 

2 4R 

5 8R 

6 


2 5R 

2 

51 :r 

6 



22R 

2 4R 

5 hr 

6 



22R 

2 

5 14R 

6 



1 7R 

2 6R 

514R 

6 



17R 

2 

51 9R 

6 



13R 

2 5R 

519R 

6 



i 3R 

2 2R 

522R 

6 



7 R 

2 8R 

522R 

6 



7R 

2 

529R 

6 



8R 

529R 

6 





44X BY 45 Y GEOMETRY - MESH AND MATERIAL SPECIFICATIONS - POISON OUT 
WITH EXPLICIT VOIDS BETWEEN DRUMS AND CORE 

4* 0.0 


1 ,08 

2.07 

3. 3 

4.2 

4.85 

5.4 

5.55 

6.2 

6. 7 

7.0172 

8.7 

9.32 

15.5894 

11. 

11.25 

11.3 

11.6712 . 

12.5 

i 3. 

13.41 

14.4 

15. 

15.6 

16.3 

17.42 

17.6 

18.25 

18.8 

19.27 

20.1 

20.98 

21.67 

22.4 

22.94 

23.7 

24.4 

24.7 

25.27 

26.16 

26.92 

27.69 

28.24 

28.6 

29.291 
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v.nO' nj m ivj r\> in & o O' m o> 



2 * 


0,0 


1,08 

2.215 

3.3 

3.45 

4.4 

4. 8 

5.54 

5.89 

6.4 

6.65 

7.2 

7.44 

7.5 

8.8 

9.07 

9. 65 

11.07 

11.84 

12.2 

13.25 

13. 88 

14.3 

14.6 

15.4 

15.7 

25.95 

16.5 

17.22 

17.6 

18.29 

19.2 

19.89 

21.41 

21.95 

22.81 

23.27 

24.43 

24.97 

25.63 

26.04 

26.62 

27.1 

27.81 

28.6 

29.291 





8 * 


30R 

1 13R 

2 

5 





27R 

1 16R 

2 

5 





24R 

119R 

2 

5 





2 1R 

1 5R 

717R 

2 

5 




2 1R 

1 2R 

72 OR 

2 

5 




18R 

1 3R 

722R 

2 

5 




18R 

1 2R 

723R 

2 

G 




13R 

1 4R 

7 

2 2R 

714R 

2 


3 

8R 

2 

5 






13R 

1 4R 

5 

2 

716R 

2 

3 

8R 

2 

13R 

12 1R 

2 

3 8R 

2 

5 



13R 

120R 

2 4R 

3 4R 

2 3R 

5 



13R 

1 2R 

2 2R 

716R 

2 4R 

3 

4R 

2 

3R 

5 







13R 

1 2R 

2 2R 

715R 

2 5R 

3 

4R 

2 


5 2R 

6 






13R 

1 

2 

717R 

2 5R 

3 

4R 

2 


5 2R 

6 






13R 

1 

718R 

2 5R 

3 4R 

2 


5 

2R 

6 







13R 

i 

717R 

2 7R 

3 2R 

2 

2R 

5 

2R 

6 







13R 

1 

71 6R 

2 8R 

3 2R 

2 


5 

3R 

6 







13R 

1 

7 1 5R 

2 9R 

3 

2 

2R 

5 

3R 

6 







13R 

T 

4 . 

715R 

2 9R 

3 

2 


5 

4R 

6 

212R 

1 3R 

713R 

2 8R 

3 

2R 

2 


5 4R 

6 






2R 

2 

71 OR 

114R 

21 OR 

-> 

3R 

5 

4R 

6 







2R 

2 

71 OR 

1 1 4R 

210R 

3 


5 

6R 

6 







3R 

2 2R 

7 8R 

U2R 

2 1 OR 

3 

2R 

2 


5 6R 

6 






4R 

2 

7 8R 

112R 

210R 

3 

2 'A 

2 


5 6R 

6 






5R 

2 4R 

7 4R 

1 12R 

2 9R 

3 


2 

3R 

5 6R 

6 






6R 

2 2R 

7 

2 4R 

112R 

2 

9R 

3 


2 3R 

5 6R 

6 





6R 

2 2R 

7 

2 4R 

11 2R 

2 

8R 

3 

2R 

2 

5 8R 

6 





UR 

2 2R 

11 1R 

2 7R 

3 3R 

2 

2R 

5 

8R 

6 







UR 

2 2R 

1 1 1R 

2 6R 

3 4R 



5 


9R 6 
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32R 

2 3R 

5 9R 

6 


32R 

2 

5 HR 

6 


3CR 

2 3R 

51 1R 

6 


3CR 

2 

513R 

6 


27R 

2 4R 

513R 

6 


27R 

2 

516R 

6 


1 OR 

414R 

2 4R 

5 16R 

6 

1 OR 

414R 

2 

519R 

6 

1 OR 

4iiR 

2 4R 

519R 

6 

7R 

4 x i R 

2 4R 

522R 

6 

7R 

41 1R 

2 

525R 

6 

5 R 

4 8R 

2 6R 

525R 

6 

5R 

4 8R 

2 2R 

529R 

6 

3R 

4 8R 

2 4R 

529R 

6 

5R 

2 7R 

532R 

6 


7R 

537R 

6 
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TABLE I. - ENERGY-GROUP STRUCTURE 


Lower energy bound a 

Number of groups 

26 

13 

4 , 

Cross- 

section set 

6.065 MeV 

1 

__ 


3.678 MeV 

2 

1 

-- 

2. 725 MeV 

3 

-- 

-- 

2.231 MeV 

4 

2 


1. 827 MeV 

5 


— 

1.353 MeV 

6 

3 

-- 

1.003 MeV 

7 

-- 

-- 

. 821 MeV 

8 

4 

1 

. 608 MeV 

9 

-- 

__ : 

.498 MeV 

10 

5 

-- 

. 369 MeV 

11 

__ 

-- 

. 302 MeV 

12 

-- 

— 

.247 MeV 

13 

-- 

— 

. 183 MeV 

14 

6 

2 

. 150 MeV 

15 

-- 

-- 

. Ill MeV 

16 

7 

-- 

67. 4 keV 

17 

-- 

-- 

40. 9 keV 

18 

8 

3 

24. 8 keV 

19 

— 

-- 

15.0 keV 

20 

9 

-- 

9. 1 keV 

21 

-- 

-- 

5. 5 keV 

22 

10 

-- 

2.0 keV 

23 

— 

-- 

748.0 eV 

24 

11 


78.9 eV 

25 

12 

-- 

.414 eV 

26 

13 

4 


a Upper energy boundary is 14.9 MeV. 


TABLE H. - ENDF/B CROSS 
SECTIONS USED IN 
GAM-H PROGRAM 


Material 

Isotope 

1006 

Lithium -7 

1009 

Boron- 10 

1140 

Carbon-12 

1013 

Oxygen -16 

1015 

Aluminum 

1111 

Molybdenum 

1126 

Tantalum -181 

1044 

Uranium -235 

1103 

Uranium -238 

1122 

Iron 





TABLE m. - MATERIAL SPECIFICATIONS BY REGION FOR RZ MODEL 


Region 

Density, a 
g/cm 3 

Region 

I 

n 

m 

IV 

V 

VI 

vn 

vin 

rx 

X 

Volume, percent 

b KelF 

b 4 c 

A1 

Mo 

Ta 

C U 

^Neoprene 

e Epoxy 

f 7 L i 3 N 
304 Stainless 
steel 

1.0 

2.44 

2.7 
10.2 
16.6 
18.9 

1.23 

3.7 
.79 

7,9 

0.467 

0. 133 










25.7 

3.62 

54.3 








4.63 
54.9 
2. 14 
7.53 

90.8 


74.9 

7.49 

4.63 
76.4 
2. 14 

42.0 

3.9 

23.0 

4.74 

54.3 

7.49 

26.3 

100.0 






23.7 

1.6 

— 









32.4 

9. 27 












77.0 












a Use of these density values with the volume percents listed to calculate atom densities will 
reproduce those used in the calculations; the density and volume percents are not absolute. 
'’Coating on fuel rods; important constituent is 1. 11-wt. % H. 
c 93. 145-wt. % 235 U, remainder assumed to be 238 U. 
d C 4 H 5 Cl. 

^ 25 ^ 34 ^ 7 ^ 2 " 

f 95. 62-wt. % 7 Li 3 N, 1.82-wt. % 7 Li 2 , 2. 56-wt. % 7 LiOH. 


TABLE IV. - EFFECT OF INPUT SPECIFICA- 
TIONS ON MULTIPLICATION FACTOR 


[Gam cross sections. ] 


S n P Z 

Number of groups 

Geometry 

RZ 

XY 

Multiplication factor 

S 4 P 0 

4 

1.041 

1.036 

S 4 P 1 

4 

1.039 

1.036 

S 8 P 1 

4 

1.036 

1.036 

S 4 P 1 

13 

1.030 

1.023 

S 4 P 1 

26 

1.029 

1.023 

Experiment 

1.011 





TABLE V. - 26-GROUP BORON-IO 


ABSORPTION CROSS SECTIONS 


[Averaged over molybdenum reflector spectrum . ] 


Group 

Cross 

-section set 

Lower energy bound a 

GAM 

GAM(ENDF/B 

1009) 

Cross section, barn 

1 

0.0822 

0.3177 

6.065 MeV 

2 

.1449 

.3914 

3. 678 MeV 

3 

.1781 

.3431 

2. 725 MeV 

4 

.2220 

.3730 

2.231 MeV 

5 

.3475 

. 5138 

1. 827 MeV 

6 

.2853 

.4112 

1.353 MeV 

7 

.1967 

.3005 

1.003 MeV 

8 

.2410 

.3421 

. 821 MeV 

9 

.3774 

.4882 

. 608 MeV 

10 

.5669 

.7144 

.498 MeV 

11 

.6476 

.8523 

. 369 MeV 

12 

.7991 

.9376 

. 302 MeV 

13 

.9014 

1.211 

. 247 MeV 

14 

1.101 

1.454 

. 183 MeV 

15 

1.358 

1.622 

. 150 MeV 

16 

1.670 

1.795 

. Ill MeV 

17 

2.209 

2. 121 

67. 4 keV 

18 

2.908 

2.660 

40. 9 keV 

19 

3.608 

3.392 

24.8 keV 

20 

4.543 

4.399 

15. 0 keV 

21 

5.819 

5.559 

9. 1 keV 

22 

7.458 

7.127 

5. 5 keV 

23 

9.872 

9. 668 

2. 0 keV 

24 

15.84 

15.87 

748. 0 eV 

25 

25.21 

25.24 

78.9 eV 

26 

75.00 

75.07 

. 414 eV 


a Upper energy boundary is 14.9 MeV. 




TABLE VI. - CALCULATED REACTIVITY WORTH 


OF CONTROL DRUMS 


Case 

Calculation 

Cross section 

Reactivity worth, 
$ 

1 

S 4 P 0 4 -group 

GAM 

9.59 

2 

4 -group 

GAM 

9.69 

3 

S 4 P 1 13 -group 

GAM 

9.56 

4 

S 4 P 1 13 -group 

GAM(ENDF) 

10.72 

5 

S 4 P 4 13 -group 

GAM except 

GAM(ENDF) 10 B 

10.45 

6 

S 4 Pj 13-group 

GAM(ENDF) 

11.61 


with explicit VOID a 



7 

S 4 P 1 26 -group 

GAM except 

GAM(ENDF) 10 B 

10.54 

Experiment (inverse kinetics) 0 

11.58 


a XY geometries 3 and 4 (appendix) used for this case; all other 
cases use' geometries 1 and 2 to obtain control swing. 
^Inverse counting method gave $11.98. 


TABLE VII. - ATOM DENSITIES FOR XY CALCULATIONS 


[Region numbers correspond to those in geometry specifications. ] 


Isotope or 
element 

Region a 

1 

(Fuel cells) 

2 

(Reflector) 

3 and 4 

(Boron carbide) 

5 

(Pressure vessel) 

Atom density, atoms/barn-cm 

H 

7 Li 

O 

N 

Al 

Mo 

Ta 

B 

C 

235y 

238y 

0.0001947 

.0129834 

.0002579 

.0042108 













b 0. 0024123 
b . 054 1190 






.0041422 


0.044116 


0.0739140 

.0184876 




.0118576 

.0008727 










a A sixth region located outside the pressure vessel is void and zero 
cross sections were used. 

b For the two cases which had an explicit void between core and con- 
trol drums (designated region 7), the reflector atom densities 
were adjusted accordingly to 0. 055586 for Mo and to 0. 0024753 for 
Al. The transverse leakage from region 7 is based on use of the 
transport cross section of region 2. 






Figurel. - Fuel element for critical experiments. Outside diameter, 2. 16 centimeters. 
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Stationary reflector^ 



Figure 2. - Cross section of critical assembly. 
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(b) XY calculation model showing region outlines with explicit void representation. 
Figure 3. - XY calculation model. Control drums in poison out position. 











1.20 1.19 1.18 1.15 1.10 .99 


.28 .11 

.29 .12 


Figure 6. - Relative power distribution in boron-carbide-controlled reactor - S^Pj 13-group XY calculation. 
Multiplication factors: k ex p° 1.011; * 1.018. 


NASA-Langley, 1972 
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